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Oscillations of Bose condensates in a one-dimensional optical superlattice 

Chou-Chun Huang and Wen-Chin Wu 

Department of Physics, National Taiwan Normal Univesity, Taipei 11650, Taiwan 

Oscillations of atomic Bose-Einstein condensates in a ID optical lattice with a two-point basis is 
investigated. In the low-frequency regime, four branches of modes are resolved, that correspond to 
the transverse in-phase and out-of-phase breathing modes, and the longitudinal acoustic and optical 
phonon modes of the condensates. Dispersions of these modes depend intimately on the values of 
two intersite Josephson tunneling strengths, Ji and J2, and the on-site repulsion U between the 
atoms. Observation of these mode dispersions is thus a direct way to access them. 
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Cold atoms in optical lattices is currently an ac- 
tive research topics in the cross-disciplined field of 
atomic molecular and optics physics and condensed- 
matter physics. From the condensed matter point of 
view, because (i) laser power can be tuned to vary the 
ratio of on-site repulsion to inter-site coupling of cold 
atoms, (ii) different geometry of laser beams can be taken 
to manipulate different dimensions and configurations of 
optical lattices, (hi) different types of atoms (fermions, 
bosons, ions, or their mixtures) can be loaded into the 
optical lattice, consequently the simulation of condensed 
matter environment is considerably easy and flexible. Re- 
cent investigations on the cold-atom system with an op- 
tical lattice have been made in various aspects, includ- 
ing the superfluid (SF)-Mott insulator quantum phase 
transition^, the band-structure phenomenal, and the 
quantum informatics^. 

In this paper, we consider a physical environment very 
common in crystalline solids. It is often in a solid that 
the system is best described by a superlattice (i.e., lat- 
tice with a basis) . This occurs when there are more than 
one type of atoms (or sites) within a unit cell. Corre- 
spondingly in the cold-atom system, a ID optical lattice 
with n-point basis can be set up as follows. One can 
shine two laser beams coherently towards the same direc- 
tion (say, z-axis), with one having n times of frequency 
(and thus n times of wave vector) of the other. Along 
with their reflected beams, the resulting electric field is 
E(z,t) =Ai cos(nkz)e inut + A 2 cos(kz)e iuJl c + c.c. and the 
dipolc energy, which is proportional to the time average 
of electric field square, is thus 



U(z) oc E(z,t) 2 ~ A\ cos\nkz) + A z 2 cos\kz). (1) 

The cross term vanishes due to time average. Conse- 
quently Bose condensates, which is prc-confined in the 
x-y plane by the magnetic trap, can be optically trapped 
and form a z-dircction ID optical lattice with an n-point 
basis. Since the ratio of amplitude A\/A% can be ad- 
justed and the factor n can be chosen, ID optical su- 
perlattice can be created for a great flexibility. Figure 
illustrates how the ID optical superlattice with a 2 and 
3-point basis is formed and their resulting potential en- 
ergies (Ai/A 2 = v5 is taken). Josephson coupling J's 
that differ due to condensates tunnelling across different 



potential barriers are also indicated. 
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FIG. 1: (a) Illustration of how a ID optical lattice with an 
n-point basis is formed by two laser beams and their reflected 
beams. Frame (b) and (c) show the resulting potential ener- 
gies (thick solid lines) as a function of lattice spacing of the 
n = 2 and n = 3 case, c denotes the lattice constant and J's 
denote the Josephson tunnelling couplings. In frame (b), d 
denotes the spacing between any potential maximum and its 
nearby minimum. 

Since the achievement of atomic BEC in 1995, dynam- 
ics of single Bose condensate had been under detailed 
investigation both theoretically and experimentally. 
When condensates are distributed in an optical lattice, 
in addition to collective modes, phonons can also propa- 
gate among the lattice. How the collective and phonon 
modes are coupled in such a condensed system thus 
has attracted a large number of studies^SiSiiSiiiii^. In 
fact, condensate excitations in a 3D optical lattice have 
been measured recently in the SF regime^. Among dif- 
ferent theoretical approaches, Martikainen and Stoo£i£ 
have recently applied the variational approach to study 
the transverse breathing mode and longitudinal acoustic 
phonon in a simple ID Bose-condcnscd optical lattice. 
The advantage of this approach is that it enables nat- 
urally the coupling between the breathing and phonon 
modes, and at the same time, allows for analytical re- 
sults. In the present context, we shall generalize their 
approach to study how the breathing and phonon modes 



2 



of Bosc condensates are coupled and behave in a ID op- 
tical superlatticc. 

At sufficiently low temperatures, dynamics of Bose 
condensates is governed by the time-dependent Gross- 
Pitaveskii (GP) energy functional 



E = / d?r\-— **V^ + [T4 rap (r) 



+ 



1*1 



(2) 



where * = ^(r,t) is the time-dependent wave func- 
tion of the condensate, p is the chemical potential, and 
Up, = Airft^a/m with a the s-wave scattering length of 
the two-body interaction. In the present case of inter- 
est, the external trap potential has two contributions, 
Vtrap(r) = V m (r)+V (r), where V^r) = ±m[uj 2 p 2 +u; 2 z 2 } 
with p 2 = x 2 + y 2 is due to the magnetic trap, while 
V {r) = V (z) = 14 cos 2 + y 2 cos 2 (^) is due to 

the optical trap. When u> p 3> lu z for the magnetic trap, 
optical potential then results a ID optical superlatticc 
(along z axis) with a two-point basis. Here c corresponds 
to the lattice constant which is related to the wavelength 
(A) or the wavevector (k) of the laser beam with fre- 
quency to by c = A/2 = ir/k (sec Fig.^l. It is convenient 
to work with dimcnsionlcss quantities and rescale the en- 



ergy, time, and length as E/(hu 



E, tujp — > t, and 



r [ip = y/K/mu>p). Eq. (J2J thus becomes 



d 3 r 



1 



2 
*| 2 



(x 2 +2/ 2 ) 



(3) 



where the dimensionless coupling V (z) = V (z) / '(Hlj p ) , 
Uo = 4ira/£ p , and p = fj,/(fko p ). The frequency u> z is 
completely dropped in ®. 

For longitudinal phonon modes to be considered un- 
der the harmonic approximation, tight-binding (TB) limit 
is applied. This means that anharmonic effect is neg- 
ligibly small. This also means that optical potential 
barriers need to be large, so condensate wave functions 
are strongly localized in z-direction around the potential 
minimums. In this limit, ^(r,t) can be taken to be the 
form 

*(r,t) = ^{ Wl [z - tc- d]^{x,y,l;t) 
t 

+ w 2 [z-(l + l)c + d]$2(x,y,£;t)}, (4) 

where £ sums over all unit cells and d is the spacing 
between any potential maximum and its nearby mini- 
mum (see Fig. QJ. Pertaining to site 1 and site 2, W\ 
and W2 are the strongly localized Wannier functions de- 
pending on z only, while $i and $2 are the ones associ- 
ated with x, y coordinates which in turn contain the peri- 
odic factor ~ cxp(i£kc) for any propagating wave (along 



the lattice direction) of wave vector k. It is noted that 
w\{z) = w 2 (—z) for the reflection symmetry. Moreover, 
for simplicity, fluctuations (time dependence) are consid- 
ered to be through &i(x,y,£;t) only. 
Substitution of J3J reduces © into 

£ = EE/ d 2 p(-^-f*(£)V 2 ^(£) 

l 1=1,2"' ^ 



+ 



E 
E 



d 2 p[J 1 <t>* 1 (£)$2(£-l) 



d 2 p{j2<s>* 1 (e)<f>2(e) + c.c. 



where &i(£) = $i(x,y,£;t) and 

2dz 2 



Jl = — J dz w*[(z — d) 
J2 = — / dz w\(z — d) 



+ V (z) 



Id 2 

-2W + V °M 



U = Uo dz |wi(2)| 4 = U dz \w 2 (z 



(5) 

w 2 (z + d) 
w 2 (z ~ c + d) 
(6) 



In ijSJ and (jBJ, only nearest-neighbor intersite couplings 
are considered. J\ and J2 are thus the two differ- 
ent nearest-neighbor Josephson couplings responsible for 
condensate tunneling. In principle, w\,W2 can be solved 
numerically, which in turn solve Ji, J2, and U. Typical 
value of Ji or J 2 is less than 0.1 in the TB limit. 

In the following, as mentioned before, we will focus on 
the breathing and phonon modes. Similar to Martikaincn 
and Stoo£i£, a Gaussian ansatz 



x exp 



B [l + eUt)+ie'Ut)}(x 2 +y 2 ) 



is assumed for the wave functions, where N represents the 
average number of atoms per site and Bq represents the 
condensate size in the equilibrium state. In fact, when 
N is fixed, the value of Bo can be calculated through 
minimizing the GP energy in . This was done in Ref 
where B Q = \ j\f\ + W with U' = 4NU/ir is given. No 
cross (xy) term is considered in (jJJ) because breathing 
mode is isotropic in the x-y plane. For each site (£, i), 
dimensionless dynamical variable e^j(i), 2 - (i) , <5^(i), 
and V£ t i(t) corresponds respectively to fluctuations of the 
local amplitude, the local phase, the number of atoms, and 
the global phase of the condensate. 

One crucial aspect on (J7J is however that it enables 
naturally the coupling between the transverse breath- 
ing and longitudinal phonon modest. In a single Bosc- 
condensed system with repulsive interaction, the fluc- 
tuation of condensed atom number is coupled to the 



3 



fluctuation of condensate size. When condensates are 
distributed in an optical lattice, the fluctuating degrees 
of freedom are complicated by the intersite coupling of 
Josephson tunnelling. For the present optical superlat- 
tice, the case is further complicated by the possible out- 
of-phase motion in addition to usual in-phase one. How 
important are the various couplings and how the conden- 
sates arc fluctuating in such a system are thus the subject 
of the following studies. 

By variational approach, one starts from the La- 
grangian of the system, L = T — E. Here T = 
J dr i^{^*d^/dt~^d^*/dt) and E is given by © [and 
hence JJJ]. Applying Q in T and E and expanding to 
second order in dynamical variables, one obtains 



-=yy 

i i=l,2 



(i + <M 



f (1 - 4,<) 



and 



-=yy 

l i=l,2 



(i + <M 



2k ( 2 + ( <* )2 ) 



(8) 



(9) 



+ : y« l ) 2 + ^o(l + 4 l )^-M 



E, 



where the Josephson-tunncling term 



+4(1/41 - ut-1,2) 2 + (&e,i - fa-1,2) 2 

-4(^_i - Vt-X,2){^Ll - 4-1,2) - 4(^,i + ^-1,2) 



2 \ -« 



4- 

1 

+4(Z^! - 

-4(1/41 - 



[(4i-4 2 ) 2 + 2(e^i-e« f 



v i,2) 

^,2)(e",i 



($1,1 ~ t>4 2 ) 2 



i)-^/,i+y. (10) 



It is noted in @ and (|10|l that when chemical potential 
\x is chosen to be p, = i/2Bo — Bo/2 — Ji — J2 (as given 
by the stability condition of dE/d8i,i\ vax —o = 0), linear 
in Su terms vanish. In this case, only second order terms 
are left with Ej in l|10|) and it justifies the smallness of 
Ej in accordance with TB. 

With (|5|)- (|10[) . one can then derive the Euler-Lagrange 
equations of motion for all eight dynamical variables. 
They can be linearized and written as 

e^i = 2B e",i + Ji( e e,i ~ ^-1,2) + -M^i _ 4',2)> 



4'+i,i) + ^(4',2 - f 
2U'B{)8(s — J\(e t 1 — 4-1, 2) 



-2Boe i2 + Ji(4',2 

-^2(641 - ££,2)1 

2 , , , 

~"5~4.2 _ %U Bo6t 2 - <M4, 2 
Bo 



4+i,i) 



-«M4,2 - 4,i): 
^41 = Ji[2(i/4i - 1^-1,2) - (641 - 4-1,2)] 

+J a [2(i*, 1 - 1/42) -(e'ir-e'iz)}, 
£42 = ^1(2(1/42 — 14+1,1) _ (4',2 _ 4+l,l)] 

+J2[2(r/42-^,i)-(4',2-4'i)]> 
£41 = -(l + 3C/')5o4,i-3S C/'(54i 

•h 
2 

J2 

vi.2 = -(l + 3J7')-Bo4,2- 3BoU'St,2 
•h 
2 

J_2 

2 



[(^,1 - ^-1,2) + (4,i - 4-1,2)] 
[(£41 - £42) 

+ 3i 

7r[(<5£,2 - ^+1,1) + (4,2 - 4+i,i] 

[(^,2-M + (4, 2 -4,i)]- 



(ii) 
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FIG. 2: Dispersions of four-branch modes resulting from flj . 
The parameters Ji = 0.09 and J2 = 0.1 are the same for 
all three frames, while U' = 1, 10, and 100 respectively for 
frame (a), (b), and (c). Here IB, OB, AP, and OP denote for 
in-phase breathing, out-of-phase breathing, acoustic phonon, 
and optical phonon modes. See text for more description. 

Note that the above eight linear first-order differential 
equations can be cast into four second-order ones in terms 
of any two sets of variables e' ti {t), e"^), 8i^{t). and 
vt,i(f) (i = 1>2) only. (That is in the lowest-order limit, 
the other two sets will share the same excitation spec- 
tra.) Considering the coupled equations between e^it) 
and #4i(i), for example, and searching for solutions of 
the type: e^(i) = e '. e K^c-^t) and 5 u {t) = 5 l e i ^ kc - UJt \ 
a 4 x 4 dynamical matrix (not shown) will be attained, 
which is then diagonalizcd to obtain four branches of 
modes. In Fig. [21 dispersions of the four branches of 
modes are shown for three cases: with same J\ = 0.09 
and J 2 = 0.1, but different U' = 1, 10, and 100. In 
the study of Bose-Hubbard model on the SF-Mott insu- 
lator quantum phase transition in 3D, it has been estab- 
lished that U/J = zx 5.8 (with z the number of nearest- 
neighbor sites) is a critical point for the case of one atom 
per site on average [N = 1). Below (above) it the system 
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is in the SF (insulating) phased. Here we are interested 
in collective excitation of condensates with average atom 
number per site being N ~ 10 — 100. Thus the system 
being in the SF phase is justified when U'/J < 100 (re- 
call V ~ NU). Consequently the U' = 1 and 10 cases 
in Fig. El are considered to be in the SF phase, while the 
U' = 100 case may be close to the phase boundary or in 
the insulating phase. 

The four branches of modes are labelled by in-phasc 
breathing (IB), out-of-phase breathing (OB), acoustic 
phonon (AP), and optical phonon (OP) modes respec- 
tively (see Fig. EJ). Strictly speaking, these descriptions 
are exact only when phonon modes are completely de- 
coupled from breathing modes and at the same time, the 
propagating wave comprised of e' e 1 (t), for instance, has 
a phase angle ir shift (out-of-phase) by the the propagat- 
ing wave comprised of e' e 2 (t). The coupling effect, which 
is proportional to the magnitude of tunnelling strength 
J, is indeed small in the TB limit. One should be noted 
that the coupling effect is also accompanied by the factor 
of BqU' (recall that Bq corresponds to the size of con- 
densate), to which it becomes more and more important 
in the large V limit. 

At the long- wavelength limit (kc « 1), dispersions of 
the four branches of modes can be analytically solved and 
uniformly written as 



2 , 2 ; 2 

= u 0i + Uf k 



(12) 



where, at the TB limit (keeping to J linear order), one 
obtains 



/ B 2 U' 2 \ 
4 + 8B (,h + J 2 )[1 + U' + -2__ J i = OB 



Oi 



8B Q U'(.h + J 2 ) ( 1 




BlU' 



i = IB 
i = OP 
i = AP, 



2B Q U'c 2 



J1J2 
J 1 + o h 



(14) 



(13) 



As shown in (|12fl -114 [) . the slopes (itj) and the k = in- 
tercept (ujQi) of the four branches of modes depend cru- 
cially on the values of Ji, J%, and U' . Therefore ob- 
servation of these modes will be a direct way to access 
them. Some features worth noting are as follows. For 
IB mode, the k = frequency (u = 2lu p ) is robust re- 
gardless of the values of J\ and J2 . It is a clear result to 
the uniform (single-condensate) limit. For AP mode, on 
the other hand, the phonon velocity at small k behaves 
like u ~ \/BoU T Jc, an expected result for the current 
system of a bulk modulus ~ BqU'J. In the large U' case 
[see Fig. Hfc)], IB and OP modes start to hybridize, to 
which a small gap (not visible in the scale potted) devel- 
ops where the two branches cross. This feature may be a 
signature when the system undergoes a quantum phase 
transition from SF into insulating phase. 

In summary, oscillations of Bose condensates in a ID 
optical lattice with a two-point basis is investigated. Fo- 
cuses are made to the transverse breathing modes and 
longitudinal phonon modes. Pertaining to condensates 
in site 1 and site 2, in-phase and out-of-phase modes 
are obtained. The mode dispersion relations depend cru- 
cially on the sizes of the two intersite Josephson tunneling 
strengths J\ and J2 and the on-site repulsion U. For the 
phonons, there arc optical as well as acoustic modes, in 
a close resemblance to those in a ID crystalline solid. 
While direct observation of these in-phase and out-of- 
phase oscillations of condensates may be resolution lim- 
ited by current instruments, dispersion relations of these 
modes should be probable by inelastic scattering mea- 
surements. 
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2Bo(l + U')c 2 
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